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A calculation of the deuteron polarization observables Ay, Ayy, Axx, A^z and the 

■ 

' differential cross-section for elastic nucleon-deuteron scattering at incident deuteron 

O '■ 

' energies 270 and 880 MeV in lab is presented. A comparison of the calculations 

> ! 

O , with two different deuteron wave-functions derived from the Bonn-CD A^A^-potential 

o 

model and the dressed bag quark model is carried out. A model-independent ap- 



(N 



proach, based on an optical potential framework, is used in which a nucleon-nucleon 
T-matrix is assumed to be local and taken on the energy shell, but still depends on 
the internal nucleon momentum in a deuteron. 



> ■ I. INTRODUCTION 

o 
o 

O ' and theoreticians as one of the clue tasks in few-nucleon physics. For three decades it has 



The reaction of elastic nucleon-deuteron scattering is considered both by experimentalists 



, been served as a hope to obtain more information about the intermediate- and short-range 
O ■ 

NN interaction and as a probe of the deuteron structure at small distances (for a review, 
^ . see Ref. During the last decade it is also studied with the purpose of testing of various 
_cd \ three- nucleon forces (3NF) and particularly their spin-dependence |2|, yj . It is also of interest 
as a basic reaction to establish a polarimetry for vector-tensor mixed polarized beams. 

To describe Nd elastic scattering below the pion production threshold different techniques 
have been applied Mia . The momentum space Faddeev equations can now be solved 
with high accuracy for the most modern two- and three-nucleon forces. The 3NFs in such 
a calculations are of Fujita-Miyazawa [sl, [ol or Tucson-Melbourne [l^ type. It was found 
from these calculations that the differential cross section and the polarization observables are 
essentially insensitive to the choice of the two-nucleon interaction provided it is in agreement 
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with A^A^ elastic scattering data. 

The next step is to go to the higher momentum transfer region and to explore this reac- 
tion above the A^A^ inelastic threshold. But calculations at incident proton energies greater 
than 200 MeV in lab encounter with some nontrivial difficulties. The ffist problem is that 
up to now there is no reliable quantitative model for the A^A^-interaction above the inelastic 
threshold. All existing models that pretend to description of two-nucleon scattering up to 



1 GeV do it only in semi- quantitative manner [11|, [1^. The second problem is that it is no 
longer conventional three-body problem. Nucleon isobar and meson degrees of freedom start 
to play a significant role and disguise the effects from the short-range two-nucleon interac- 
tion ^]. And last but not least, there are purely computational difficulties in performing 
relativistic Faddeev calculations at higher energies Hence, the existing theoretical 

frameworks of elastic Nd scattering at intermediate energies concern mainly to the very 
backward scattering angles where short-range behavior of a deuteron wave function (DWF) 
is, as expected, the most evident. The common framework at these energies is a multiple 
scattering formalism, where the two-body Nd amplitude is expanded in series of a A^A^ T- 
matrix. In the absence of a reliable two-nucleon potential in the GeV region, the T-matrix 
is usually taken either on-energy-shell |16] or extrapolated off-shell in a model independent 



171 ] . So, the DWF contains the only information about the A^A^ interaction in these 



fashion 
models. 

Now, there is a set of high quality nucleon-nucleon potentials, based mainly on a meson 
exchange picture, that describe deuteron properties and two-nucleon elastic scattering data 

n n □ 

perfectly below the pion production threshold p^, 119|, 120|. With the inclusion of a 3NF, 
they also provide a good description for the 3N binding energies and the nucleon-deuteron 
differential cross section. The general feature of these models is an uniform depletion of the 
DWF at small internucleon distances. 

In the last two decades, numerous attempts have been made to describe the intermediate 



and short-range A^A^ interaction by methods of QCD 2l|, |22l, l23|, l2J, |25|] . The advantages 



of the quark models are that only a few physically meaningful parameters can be used to 
describe processes involving hadrons and a possibility to take the underlying symmetries of 
QCD directly into account. Most of these attempts focused only on elastic A^A^ scattering 
at rather low energies and/or on low partial waves and are not suitable as input in few- 
body calculations. However, some of the features of a few-nucleon system, derived within 
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the frameworks, can be further explored and tested. Especially it concerns the short-range 



behavior of the deuteron and 3N bound states. One of the successful models 
six-quark (6g) symmetries is a dressed bag model (DBM) proposed recently 25 



3ased on 



26|. The 



essence of this model lies in different dynamics of the quark configurations |s^p^[42]2;L = 0, 2), 
|s^p^[33]a;L = 1,3) and the most symmetric ones Is^fOj^^L = 0), |s^p"'^[51]a;Iv = 1). Whereas 
the first two configurations have a cluster structure and are projected mainly on the A^A^ 
channel, the third and fourth configurations have the structure of a quark bag with a large 
weight of the AA and CC (hidden-color) states. As a result, the deuteron and 3N wave 
functions have a short-range node, which result from the orthogonality of the cluster wave 
function and the wave function of the 6g compound state. These 6o symmetry arguments 



are also the foundations for the Moscow A^A^-potential model [27|, l28j. In the framework of 
the DBM, it turned out to be possible to fit very reasonably the A^A^ phase shifts in ^5*0 and 
^5*1 Di channels up to 1 GeV and the deuteron static properties as well. As to the 3N 
systems, this model could explain quantitatively all static properties of ^He and ground 
states, including a precise parameter-free description of the Coulomb displacement energy 
of ^He — and all the charge distributions in these nuclei 29|, l30l |. 

In this work, elastic Nd scattering is considered as a possible discriminative tool between 
calculations with two different kinds of a DWF. The first one is derived from the Bonn-CD 
A^A^ potential 20| and diminishes uniformly with approaching the internucleon distance to 



zero. The second one is a result of the DBM [25|] and develops a nodal behavior. The 
framework is based on an optical potential formalism. A multiple scattering expansion 
is used to derive the optical potential that comprises of the one-nucleon-exchange (ONE) 
mechanism, the single- and double-scattering terms. As shown by Faddeev calculations, 
at energies below about 200 MeV rescattering of higher order is very important, however 
around 300 MeV the first two terms in the expansion are sufficient to describe the total Nd 
cross section 311] • For the A^A^ input, a model-independent approach is used in which the 
nucleon-nucleon T-matrix is assumed to be local and taken on-shell, but in contrast to the 
common impulse approximation it still depends on the internal nucleon momentum. Thus 
the two-nucleon amplitude cannot be factorized out of an integral and some kinematical 
off-shell effects are implicitly taken into account. 

The deuteron vector and tensor analyzing powers and the differential cross section are 
calculated for two deuteron kinetic energies Ed in lab - 270 MeV and 880 MeV. At Ed = 270 
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MeV a comparison with precise data measured at RIKEN 32] is performed to validate 
the model. At this energy the difference between the calculations with the two kinds of a 
DWF is, as expected, not remarkable. The calculation at 880 Mev is performed in view 
of the recent experiment at JINR 3^, in which the deuteron polarization observables were 
measured in the region of the so called "cross-section minimum" {6cm = 70 — 140°) where 
both vector and tensor analyzing powers can have large values. 

The structure of the paper is as follows. In Section 2, the formalism of the optical potential 
and multiple series expansion is briefly given. The approximation to the fully off-shell A^A^ 
T-matrix is developed and the Nd scattering amplitude is defined. Section 3 is devoted to 
results of the calculation and their discussion. Conclusion summarizes the content of the 
work. 



II. THEORETICAL FRAMEWORK 



The basic equations for a three-body system in quantum physics are Faddeev equations 
which can be written in the operator form 

U = PG^^ + PTGof/, (1) 

where U = U^'^^'^^^^^j^{q',q) is an amplitude of elastic A^cZ-scattering, q (g') - initial (final) 
relative momentum in the nucleon-deuteron cm., fidjfJ'N are spin quantum numbers, T - 
A^A^ scattering T-matrix, Gq = {E ~ Hq + ie)^^ is a free propagator of the 3A^ system and 
P = P12-P23 + -P13-P23 stands for a permutation operator that takes into account the property 
of identity of three nucleons. 

One can rewrite ([T]) in the form that is more appropriate for calculations at higher energies 
where a multiple scattering expansion is justified [3J]: 

U = Kpt + VoptGdU. (2) 

This is an optical potential framework and is a nucleon-deuteron optical potential: 

Kpt = PG'o^ + PT,GoKpt. (3) 

Here is a A^A^ T-matrix without the deuteron pole term and G^ is a deuteron contribution 
in the spectral decomposition of the two-body Hamiltonian. 






(a) 



(b) 



FIG. 1: The Nd optical potential up to the second order in the multiple series expansion: 

(a) - one-nucleon exchange, (b) - single scattering, (c) - double scattering; dark-filled circles 

represent the A^A^ T-matrix and grey-filled - deuteron wave function. 



Although, Eq. ([2]) seems quite simple and represents a common two-body scattering 
equation, the derivation of the optical potential from Eq. (j3]) is as difficult as solving the 
Faddeev equations themselves. However, one can employ a multiple scattering expansion of 
Vopt, keeping a possibility to include in Vopt some additional terms originated, for example, 
due to a 3NF. At hi ghe r energies only a few terms in this expansion may be sufficient to 



describe observables 



3l|: 



Kpt = PGo' + PT,P + PT,GoPT,P + 



(4) 



In Fig. [T] a schematic representation for the optical potential (up to the second order in 
series expansion) is shown. The first diagram (Fig. [U^a)) is a usual one-nucleon exchange 
mechanism, the second graph (Fig. [T](b)) is a so called triangle diagram and represents the 
single-scattering process. And Fig. [U^c) represents the double-scattering which takes into 
account the break up channel and the double-charge exchange mechanism. 

If to define an initial state as \i) = \q; yUd/iAr)i(23), where a proton 1 is in continuum and a 
proton 2 and a neutron 3 are bound in the deuteron, then (jl]) can be written in the matrix 
form: 

^iMmGoT^^''^\^) + 1(23) (/iTi^^^^^Gof.lz) + 2i3l){f\n^''^GoT^n^)■ (5) 

Here, 

1(23)(/|T, = 1(23)(/|T5^^)+2(31)(/|T5^^) 

is an antisymmetrized proton-proton T-matrix and T(f means that the interaction is oc- 
curred between the particles i and j whereas the particle k is a spectator. 
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The one-nucleon exchange contribution is written in a usual way: 
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q^ + Ml + \lAq' cos^ ^ + Mf, - + Mj 



^Is [q^ \q'^ (q' + \ 



Here, a comphcated notation related to a summation over spin and orbital quantum numbers 
is omitted, is a wave function of the deuteron composed of nucleons % and j, 6' is a 
scattering angle and 

,^ _f EdMl{Ed + 2Md) V'^ 
^ ^ \{MN + Mdf + 2MNEj ' 
where is a kinetic energy of the deuteron in lab. 

To evaluate a single scattering diagram, one must implement an integration over the 
internal momentum of nucleons in the deuteron. To do this, a knowledge about the fully 
off-shell behavior of the A^A^ T-matrix is required: 

i(23)(/l^cK) = / ^,^1, {p+ T, {q\p- \q' + \q, q,p- + \q') ^23 (p- , 

(6) 

here k = q — q' is a transferred momentum. 

However, the absence of a high quality A^A^ interaction model above the inelastic thresh- 
old at the present moment forces anyone to make some approximate evaluation of the inte- 
gral, proceeding from the assumption of either off-shell behavior of somehow parameterized 
T-matrix or taking only its on-shell value. The common approach of twenty years old calcu- 



lations is an optimal impulse approximation 
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361] in which the internal momentum in the 



T-matrix is put to zero in ([6]) that permits the T-matrix to be shifted outside the integral 
and the remaining integration produces the deuteron form factor. In this approximation the 
leading order corrections due to Fermi motion is vanished provided that the A^A^ amplitude 
is spin independent and local. Thus the variation of the A^A^ T-matrix with momentum 
p is compensated to leading order in p/M^ by appropriate choice of the energy parameter 
upon which the T-matrix depends. The condition on the energy is that the A^A^ T-matrix 
be on-shell when evaluated at p = in Eq. ([6]). However, this approximation misses some 
vital momentum dependency in the T-matrix and is suitable only for scattering on a heavy 
nucleus where the recoil is not significant. 

In this work, the T-matrix is put on-shell in a way that minimizes the off-shell corrections. 
Firstly, the A^A^ T-matrix should be transferred to the two-nucleon cm. frame, where it is 
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usually defined. By means of the Lorenz transformation /^(^(Pcm), one has: 

Tc [q'J- \q' + \q; q,p- + \q'^ = tv{C)Tr (Q', Q) A-\C), (7) 

where Q' = {C-'[q'] - [p - |g' + \q]) /2 and Q = {C-'[q\ - [p- fg + \q']) /2 
are final and initial two-nucleon relative momenta in A^A^ cm. and A is a transition operator 
that includes the boost and the Wigner spin rotations of the T-matrix; Pcm = P + \q' + \q 
- cm. momentum of the two-nucleon system. 

If the T-matrix is not strongly energy-dependent and essentially local, then it depends 
on two variables - transferred momentum {Q' — QY for the direct A^A^-interaction and 
{Q' + QY in case of the exchange mechanism. Then one can introduce two new variables 
Q and 9 that correspond to the on-shell relative momentum and scattering angle and define 
them as follows: 

|g'-g| =2Qsin(^/2), |g' + g| = 2gcos(^/2). (8) 

Thus, the T-matrix T^^™ is taken to be on-mass shell and calculated at the effective two- 
nucleon energy in lab 

Thus, although the A^A^ T-matrix is evaluated on-shell, it still contains some off-shell infor- 
mation, since the effective energy -Ee// and the scattering angle 6 depends on the off-shell 
momenta Q' and Q. Particularly, the effective energy depends on the internal nucleon mo- 
mentum p, and therefore the T-matrix cannot be driven outside the integral in Eq. (I6l). 
Moreover, some off-shell dependency is present on a relativistic level in the operator A when 
the transformation of the T-matrix from the A^A^ cm. frame to the nucleon-deuteron cm. 
frame is performed according to Eq. ([7]). In the relativistic case, this approximation corre- 
sponds to the situation when the T-matrix depends only on the two kinematical invariants - 
t and u, and its evaluation is performed via imposing the on-shell condition on the squared 
total energy s = 4M^ — u — t. 

For the double scattering term one should evaluate a six- dimensional integral over the 
internal momentum in the deuteron and the intermediate momentum of the scattered nu- 
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cleon: 

i(23)(/|r2G'oTi|i) = 

X Ti (•?" + i(4 + 9'). -P - i(9' + 9"); 9. -P - i(9-- 9"')) l-zs (p - y) , (10) 

The A^A^ T-matrices in fllOp are approximated in the same manner as for single-scattering 
according to Eqs. (I7l)-(l9]). The Green function Go is taken in a nonrelativistic form to 
simphfy the evaluation of the pole part of the integral: 

1 ^ ^ 0^ _ {q"+{q + q')/'2f _ {p-{q + q")/2Y _ ip+iq' + q")/2? 
° 2Mn 2Md 2Mn 2Mn 2Mn 

Assuming here Md = 2Mn, one can rewrite it as follows 



r.-i _ jq' - gT)/2 - R 

~ 2M^ 



(12) 



where 



R' = + 2p' + lq"iq + q') + p(g ' - q). (13) 



Then, redefine the variables of integration and instead of the variables \q"\ and |p| introduce 
two new variables R and a as 

=Y --Rcosh/?cosa, (14) 

|p| ='^^i?cosh/3sinQ;. (15) 

The expression for (3 can be found by substituting these definitions in Eq. f[T51) . Thus, the 
evaluation of the pole part of the integral ( fTOl) is straightforward. Here, the both the pole and 
the principal parts of integration in the double-scattering term is taken into account. The 
evaluation of the integrals is performed by means of the Monte-Carlo simulations, dividing 
the integration domain on several parts to minimize numerical errors. 

The scattering equation ([2]), which is here a relativistic Lippmann-Schwinger equation, is 
solved in helicity basis employing a i^-matrix approximation, i.e. only the pole part of the 
two-body propagator Gd is remained thus all terms in the equation contain only on-shell 
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information about the optical potential and the scattering amplitude: 
{X'MU'{q',q)\XAN) = {X'MVit{Q',Q)\^^^^^ 



2^ f dq" 



j ^^^2((^'MVitiQ',qlU'{q'',q)\XAN)A{q'\q^^ 



where the kinematical factor is 



A{q ,q)=q 



E?, + El, 



and Eg = a/^^ + + ^/q'^ + Mj is a total energy. 

Then for the Nd scattering amplitude in spin space one has (the incident particle is going 
along the z-axis): 

X di,^_,,^^^^^^^ie){X'MU'iq,q)\ - fi.iXN), (17) 

where J is a total angular momentum. 



III. RESULTS AND DISCUSSION 

The calculation of the deuteron polarization observables Ay, A^x, Ayy, A^z and the dif- 
ferential cross section is performed with two different deuteron wave functions. The first 



one is a wave function derived in the meson-exchange Bonn-CD model [20|. The general 
trait of wave functions of this kind, derived from the most modern A^A^ potentials, is their 
uniform behavior at small distances. The other possible choice is a wave function with a 
nodal behavior. This node corresponds to the so-called forbidden state in a A^A^ system as 
a consequence of the six-quark dynamics and the fact that the mostly symmetric six-quark 



state \s^) has a small A^A^ component 37|. As one of the representatives of such a wave 
functions with nodal behavior can serve a wave function of the DBM [25] , which has a node 
in the ^Si wave at a distance of ~ 0.6 fm. 

Although the A^A^ T-matrices in the Nd amplitude are taken on-shell, they still depend 
on the off-shell momenta, particularly on the internal nucleon momentum in the deuteron. 
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FIG. 2: Deuteron vector and tensor polarization observables at the energy = 270 MeV in lab. 
The solid and dashed curves are calculations with the DWF derived in the Bonn-CD model and 
the DBM respectively. The experimental data are taken from Ref. 32 1. 



This dependency is hidden in a value of the effective on-shell energy and A^A^ scattering 
angle. The integration on the internal momentum means that the knowledge of the T- 
matrix at a large energy interval is required. The A^A^ T-matrices are calculated using the 



recent partial wave analysis SP07 38|] which extends to 3 GeV for pp scattering and 1.3 GeV 
for np scattering. All partial waves up to the total angular momentum Jnn = are taken. 

In Figs. [2] and [3]^a) the calculations of the observables at the deuteron energy Ed = 270 
MeV are shown for the two different deuteron wave functions. The experimental data are 
taken from Ref. [s^. The convergent results in the sum (fT71) are obtained at J = 25/2. The 
full curve is a calculation with the Bonn-CD DWF, and the dashed curve is a calculation 
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FIG. 3: The differential cross section at energies = 270 and 880 MeV. The definition of the 
curves is the same as in Fig. [2j 

with the DWF from the dressed-bag model. As one can see, the two calculations do not 
differ significantly from each other. The difference between them is of the same order as 
a disagreement with the experimental data and is seemed to be caused by approximate 
treatment of the off-shell effects in the A^A^ T-matrix. Whereas for the Bonn-CD calculation 
the assumption of locality of the A^A^ potential and amplitude may be a good approximation, 
the A^A^ potential in the DBM is highly nonlocal and energy dependent, thus some off-shell 
effects from the T-matrix of the DBM may cancellate the effects from the nodal behavior 
of the DWF. As for the differential cross-section, the calculations practically coincide with 
that derived from a solution of the Faddeev equations without a 3NF. The lack of the cross- 
section at intermediate scattering angles is a common feature of such a calculations. Thus 
the higher rescattering terms in the Nd optical potential and dynamical off-shell effects are 
not significant for the cross section at this energy. The deuteron polarization observables 
are also in a good agreement with the experiment. It should be noted here, that even the 
Faddeev calculations that include a 3NF are only partially successful in a description of the 
polarization data. 

The convergence of the calculations for forthcoming experimental data at E^i = 880 MeV 
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FIG. 4: The same as in Fig. [H but for = 880 MeV. 



331] is achieved at the total angular momentum J = 39/2. This is a quite large value, hence 
the Faddeev calculations are very difficult to solve at higher energies in the partial wave basis. 
As can be seen from Fig. HI the differences between the curves derived from the two DWFs 
are remarkable. Only for small scattering angles 9 < 50°, where the transferred momentum 
is not large, the two calculations give almost the same results. It is interesting to note, that 
at very backward scattering the differential cross-section is mostly insensitive to the kind 
of a DFW and the maximum difference becomes apparent in the cross-section minimum 
region (see Fig. [3](b)). However, to make quantitative calculations at these energies the 
3NF contribution due to the excitation of the A-isobar must be taken into account. This 
contribution is the most likely mechanism to render the cross-section fall-off at backward 
angles [3]. Furthermore in the dressed-bag model, 3.6% to the DWF contributes from 
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the 6g-bag which at this scattering energies can take a large transferred momentum. So, 
the 3NF, originated from the scattering of a nucleon on this quark bag, can also provide a 
significant contribution at backward scattering angles. Anyway, such a 3NF provides a large 
amount of the and '^He bound energies [29|. 



IV. CONCLUSIONS 

A calculation of the deuteron polarization observables Ay, Ayy, A^x and A^z and the 
differential cross section in an optical potential formalism for elastic nucleon-deuteron scat- 
tering at incident deuteron energies = 270 and 880 MeV was presented. Under the 
investigation was the calculations with two different deuteron wave functions derived from 
the Bonn-CD A^A^-potential model and the QCD-motivated dressed bag model. For the 
A^A^ input, the model independent approach was employed, in which the nucleon-nucleon 
T-matrix was taken to be on-shell in a way that some kinematical off-shell effects were in- 
corporated in the definitions of the effective A^A^ scattering angle and energy. It was found 
that the differential cross section is not affected by the higher rescattering contributions and 
the off-shell effects in the A^A^ amplitude have a minor influence on this observable. At the 
energy Ed = 880 MeV, large differences in the observables calculated from the two DWFs 
were observed in the cross-section minimum region both for the analyzing powers and the 
differential cross section. However at the backward scattering angles, the differential cross 
section was shown to be mostly insensitive to the short-range deuteron structure. 

The author thanks Dr. V.P. Ladygin for helpful discussions and Prof. V.I. Kukulin 
for providing a code with the parameterized deuteron wave function of the dressed bag 
model. This work is partly supported by the Russian Foundation for Basic Research, grant 
07-02-00102a. 
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